Graphene Quantum Dots: Beyond a Dirac Billiard 
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We present realistic simulations of quantum confinement effects in ballistic graphene quantum 
dots with linear dimensions of fO to 40 nm. We determine wavefunctions and energy level statis- 
tics in the presence of disorder resulting from edge roughness, charge impurities, or short-ranged 
scatterers. Marked deviations from a simple Dirac billiard for massless fermions are found. We find 
a remarkably stable dependence of the nearest-neighbor level spacing on edge roughness suggest- 
ing that the roughness of fabricated devices can be potentially characterized by the distribution of 
measured Coulomb blockade peaks. 
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Graphene [l|, [H, the first true two-dimensional (2D) 
solid, is attracting considerable attention, mostly due 
to unique dynamics of electrons near the Fermi energy 
which closely mimics that of a massless Dirac Hamilto- 
nian [see Figfjja)]. Moreover, the double cone structure 
near the K and K' points of the sublattices in recip- 
rocal space gives rise to a near "pseudospin" degener- 
acy, suggesting an analog of Dirac four spinors. Envi- 
sioned applications range from high-mobility nanoelec- 
tronics 0] , spin-qubits in graphene quantum dots 0] and 
the creation of "neutrino" billiards 0, Dirac (in- 
cluding neutrino) billiards receive growing interest as a 
complement to classical and quantum (Schrodinger) bil- 
liards, which have taken central stage in studies eluci- 
dating the quantum-to-classical crossover in both reg- 
ular and chaotic devices. Additionally, spin coherence 
times in graphene are expected to be very long due to 
weak spin-orbit and hyperfine couplings 0, Q making 
graphene quantum dots promising for future spin based 
quantum computation [4|. However, confining electrons 
in graphene is a challenge, due to the gap-less electronic 
structure and the Klein tunneling paradox [j| [l(| . This 
difficulty has recently been overcome by structuring 2D 
graphene and quantum mechanical confinement effects 
have been observed in nanoribbons [U, [l2|, interfer- 
ence devices [l4[ , single electron transistors [lfj. flij and 
graphene quantum billiards [(|. 

In the following we present a realistic simulation for 
the single particle spectrum of graphene quantum dots 
(i.e. billiards) by explicitly considering rough edges and 
disorder. This work was motivated by recent advances 
in fabricating dots with linear dimension d ranging from 
a few hundred nm down to about 40 nm, and determin- 
ing their nearest neighbor energy level spacing distribu- 
tion We analyze dot wavefunctions, the density of 
states (DOS) and the nearest neighbor spacing distribu- 
tion (NNSD). We address the question to what extent the 
electron spectra now experimentally accessible via mea- 
surements of Coulomb blockade peaks reveal information 
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FIG. 1: (color online) (a) Graphene dispersion near K and 
K' points of the infinitely extended sheet, (b) Rectangular 
segment of a graphene flake, vertical edge armchair, horizon- 
tal edge zig-zag terminated, (c) (Approximately) rectangular 
quantum dot with rough edges. 



on the roughness and size of the graphene quantum dot. 
To put it provocatively: Can one "hear" the rugged shape 
of a drum if it is made of a graphene flake? 

We investigate graphene dots with linear dimensions 
between 10 and 40 nm, in line with currently fabricated 
devices |(|. The shape is chosen regular in the absence 
of edge roughness. One motivation of this choice was the 
remarkable result Q that a Dirac neutrino billiard, in 
sharp contrast to a Schrodinger billiard, would feature 
chaotic dynamics. We consider a ballistic dot, thereby 
neglecting inelastic scattering. This is justified as the in- 
elastic mean free path Xu found in experiment exceeds 
the linear dimension d ■■ 



yjALW/ir of the dot, d < 40nm 
400 nm [l[ . Our simulation allows for the inclu- 
sion of disorder through (i) edge roughness [see Fig.fJJc)], 
(ii) short-range disorder due to point defects in the in- 
terior, and (iii) long-range screened Coulomb distortion 
due to charge deposition. Rough edges are simulated by 
modulating the boundary of the dot by steps of height 
±5w and length ±51 randomly chosen from the interval 
[0, AW], AW < d. We refer to AW as the amplitude of 
edge roughness which varies between 0.3 nm (weak dis- 
order) and 2 nm (strong disorder). The resulting piece- 
wise straight edge features alternating zig-zag and arm- 
chair sections [see Figs. [ljb,c)]. We describe impurities 



2 



and defects in the flake by positioning either long range 
(V(r) = V e- a \ r - T °\) or short range (V(r) = 5(r - r )) 
scatterers at randomly selected lattice sites r . We use an 
impurity density rii < 1.8- impurities/carbon (10 to 
100 defects per flake), as estimated by recent work (l7| . 
The spectrum of the graphene quantum dots is deter- 
mined employing a Lanczos algorithm (l8j giving the 500 
eigenstates closest to the Fermi edge. The graphene flake 
is described by a third nearest neighbor tight-binding ap- 
proximation to correctly reproduce the graphene band- 
structure [nj. The modified C-C bond length at the 
flake boundary is accounted for by increasing nearest- 
neighbor coupling to the outmost carbon atoms by 12% 
in accordance with recent ab-initio density functional 
calculations [13]. Our ensemble averages for the DOS 
= ~ encompass typically 5000 disor- 

der realizations £. 




FIG. 2: (color online) Ensemble-averaged density of states 
(p(E)}^ of graphene quantum dots with increasing edge 
roughness, see different values for AW in the subfigures. The 
size of all devices is equal, d = 20nm (20.000 atoms), their 
width W = 16nm. Dashed lines indicate the averaged linear 
DOS for Dirac billiards [see Eq. JTJ]. The inset shows the 
K — A" splitting of 12meV. 

The linear dispersion relation of a massless Dirac par- 
ticle implies a DOS linear in e, 
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where d = ^AWL/tt is the effective diameter of a dot 
with area WL and e is measured relative to the coni- 
cal intersection [Fig. QJa)] assumed to coincide with the 
Fermi edge. The simulated DOS for the quantum dots 
display marked deviations from Eq. |lj. For weak disor- 
der pronounced size quantization peaks appear [see e.g. 
vertical arrow in Fig. EJa)]. Their positions are deter- 
mined by the smallest linear dimension of the flake [2l[ . 
We have investigated both cases W < L (W > L), and 
find the same qualitative behavior. In the following, we 
assume W < L. Note that width (W) and length (L) 
of the rectangular flake are not strictly equivalent as the 
vertical boundary features an arm-chair border while the 
horizontal forms a zig-zag border. The distance between 



the quantization peaks is AE = Hvftt/W « 1.5/W[eV], 
width W given in nm. This yields an analytic prediction 
for the energy separation of O.leV for the peaks in Fig. [21 
Weak disorder, i.e. small edge roughness can induce cou- 
pling between the cones at K and K' [Fig. [T]Ja)]. This 
manifests itself by a fine structure of size quantization by 
lifting the degeneracy [inset in Fig. [2^a)]. The quantum 
confinement peaks in the graphene dot are enhanced com- 
pared to a corresponding Schrodinger billiard of the same 
geometry in part because of both the altered dispersion 
relation and the additional degeneracy Strong disorder 
smears out size quantization patterns and the DOS be- 
gins to resemble that of a zero-mass Dirac fermion in free 
space (Eq. [1]). Only when the edge roughness can be lim- 
ited to the sub-nanometer scale, quantized conductance 
in graphene nanoribbons persists [13[ . 

Even in the limit of strong disorder, the prominent 
peak in the DOS near the Fermi edge remains unchanged 
[see horizontal arrows in Fig. Hfb)] . A direct look at the 
wavefunction [Fig. [UJa)] reveals its origin: A large num- 
ber of strongly (Anderson-) localized states at the edges 
of the graphene flake. Each eigenstate features a non- 
vanishing amplitude only at a few, not always spatially 
connected carbon atoms, with a decay length into the 
bulk of typically 0.5 nm [see left arrow in Fig. 3(a)]. We 
find that the eigenencrgies of these states are extremely 
sensitive to the site energies at the corresponding lattice 
sites. This agrees well with the experimental observa- 
tions of sharp resonance in the electron transport through 
graphene constrictions fl5j |. 

Delocalized states contributing to the size quantization 
peaks show pronounced features well beyond the simple 
picture of a confined zero-mass Dirac particle. While the 
transverse quantization resembles that of a conventional 
conductor, the interference pattern in the electron prob- 
ability density [Fig. [^b)] results from the simultaneous 
presence of multiple wavelength scales for the cone near 
the K point (unlike wavefunctions near the T point), in 
k = (fco + vpir/L, vftt/W): Parallel to armchair edges 
(i.e. in vertical direction in Fig. [3]), the wavelength is 
of the order of twice the width of the ribbon-like dot 
« 32 nm. Parallel to zigzag edges, (i.e. in horizontal di- 
rection in Fig. [3j the wave oscillations are much shorter 
with a typical wavelength of 0.7 nm [see Fig.[3]Jb)] result- 
ing from beating (frequency ratio 3:2) between lattice pe- 
riodicity a = 0.24 nm and the characteristic wavelength 
Ao = 27r/fco ~ 0.37 nm, where kg S> v-pir/L is the dis- 
tance between T and K point in reciprocal space [see 
Fig. 1(a)]. We find beating patterns with this character- 
istic length scale to be universally present in all delocal- 
ized states, even in the presence of long-range disorder. 
Only because of the sub-nanometer length scale of Ao at 
the K point is the graphene dot sensitive to edge rough- 
ness and disorder on a length scale of a few nanometers, 
in contrast to a Dirac cone at the T point. The reason 
for the AT-point (fcg,0) [as opposed to those at 60 and 
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FIG. 3: (color online) Eigenstates of graphene quantum 
dots with dot area of about 225nm 2 (13.000 carbon atoms). 
Eigenenergies are (a) 65 meV (b) 220 meV (c) 650 meV 
(d) 1020 meV respectively. The edge roughness parameter 
AW = lnm in (a-c) while m = 0.005 in (d) (White dots 
mark the center of long-range impurities). 



120 degrees, h(ko, ±V3A;o)], to appear in the eigenstate 
shown is the orientation of the zigzag (armchair) edges 
in the flake parallel (orthogonal) to the (fco,0) direction. 
However, all three directions appear for higher transverse 
quantum numbers, resulting in enhancements along the 
three zigzag-directions of the lattice [i.e. horizontal, 60 
and 120 degrees, see dashed lines in Fig. EUc)]. As a con- 
sequence, eigenstates feature a 2D hole ( "swiss-cheese" ) 
pattern emerging from the interference of plane waves ro- 
tated by 60 degrees relative to each other [see Fig. Efd)]. 

In order to delineate the influence of disorder and edge 
roughness on the energy level statistics, we have deter- 
mined the NNSD, P(AE), i.e. the probability that the 
energy difference between two adjacent energy levels is 
AE, for different amplitudes AW of roughness. Within 
the framework of quantum dynamics of Schrodinger bil- 
liards, P(AE) follows a Poisson distribution for separa- 
ble (classically regular) shapes while it should display 
a Wigncr-Dyson (or GOE) distribution for irregularly 
shaped (classically chaotic) billiards. In contrast, even 
rectangular shaped Dirac neutrino billiards have been 
shown to feature a GUE distribution because of the bro- 
ken time reversal symmetry due to chirality Af- 
ter spectral unfolding [s = p(£i), with (P(s)) = 1 and 
(sP(s)} = 1] we find for the ideal rectangular graphene 
dot (AW = 0) a near-perfect Poisson distribution [see 
Fig. HJa)] . By gradually increasing either the edge rough- 
ness or the defect concentration rii, the distribution 
smoothly evolves into a Wigner-Dyson like statistics [see 
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FIG. 4: (color online) NNSD P(s) of different billiards: Rect- 
angular graphene flake with (a) smooth edges (/3 — 0.07). 
(b) finite edge roughness AW — 0.5 nm (f3 = 0.4), (c) 
AW = 1 nm (f3 = 2.5). (d,e) Schrodinger billiards with 
same edge roughness as (a,b). The solid red curve shows fits 
to the Hasegawa distribution. Dashed (dotted) lines show a 
Wigner-Dyson (Poisson) statistic as guide to the eye. 



Fig. (Hb,c)]. Clearly, such a behavior reflects the conser- 
vation of time-reversal symmetry in graphene quantum 
dots. Among the distribution functions suggested for 
the transition regime for classically mixed phase space 
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26l |. the best fit for the disorder param- 
eters and geometries investigated was achieved for the 
two-parameter Hasegawa distribution j2?| 



P H {s;a,f3) = N 
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where p and N are determined by the normalization con 
ditions 
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271 ] . While the control param- 
eter j3 describes the transition from Poissonian (/? = 0) 
to Wigner-Dyson statistics (ft — ► 00), a is a system- 
specific constant. Indeed, we find a = 0.75 to correctly 
reproduce our numerically obtained NNSD for different 
values of both edge roughness as well as scatterers [see 
Fig. QJa-c)]. A strong edge roughness of 2nm (or im- 
purity concentration rii = 5 • 10~ 3 ) is required to reach 
the completely chaotic limit, i.e. a Wigner-Dyson NNSD 
statistics. Remarkably, for moderate values of the edge 
roughness amplitude (Ajy = 0.6 nm) a Schrodinger bil- 
liard and a graphene billiard of the same geometry dis- 
play a markedly different NNSD (Fig. g}: While the 
Schrodinger billiard has already reached the Wigner- 
Dyson limit ((3 — > 00), for the graphene the NNSD still is 
closer to the Poisson limit, pointing to the unique spec- 
tral properties of graphene. The quasi-regular dynamics 
in graphene is more stable against disorder than in corre- 
sponding Schrodinger billiards. The origin of this, at first 
glance, surprising finding is closely related to the elec- 
tronic structure of graphene at the K point [Fig. HJa)]. 
In a classical rectangular ballistic billiard with only rect- 
angular edges along the armchair or zig-zag direction 
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FIG. 5: (color online) (a) Dependence of the control parame- 
ter j3 for the transition from a Poisson to a Wigner-Dyson dis- 
tribution (Eq.[2J) on the edge roughness amplitude AW or the 
defect density m. (b) Dependence of rescaled mean-level spac- 
ing 6 = (AE ■ E/Eo)meV, where Eo was taken at lOOmeV, 
on the size of the quantum dot. Triangles represent experi- 



mental data, rescaled by AE -> {AE—E )-§- = (AE -5)^ 
to take into account the charging energy Ec and the energy 
dependence of the DOS. [f|. 



an additional constant of motion, the magnitude of the 
linear momentum |k|, exists. Such billiards are there- 
fore classically regular irrespective of the number or size 
of the edges. By contrast, due to the larger deBroglic 
wavelength of Schrodingcr billiards with wavcnumbcrs 
near the T point cannot resolve the exactly rectangularly 
shaped edges, and thus mimics chaotic dynamics. 

Of potential technological significance is the depen- 
dence of the NNSD on the disorder in graphenc billiards. 
For all three classes of disorder we consider (edge rough- 
ness, short and long range disorder) we find a linear rela- 
tion between the NNSD parameter f3 and the edge rough- 
ness amplitude, (3 m 2AW, and between (3 and the defect 
density, (3 w 0.7rij [Fig.^a)]. As (3 can be obtained with 
high precision from a fit to P(s), the edge roughness or 
defect density can be deduced from the NNSD, if the dis- 
tribution is Poisson- like in the limit AW = (rii = 0). 
We suggest that this dependence could be used to es- 
timate the disorder in experimentally realized regularly 
shaped ballistic graphene quantum dots. Numerically, we 
find /3 w 4 in a fit to data from recent experimental inves- 
tigations of a 40nm graphene billiard [6|, corresponding 
to an effective roughness AW e w 2 nm, or an effective 
defect rate n\ ~ 5.5 • 10 -3 (see black triangles in Fig. [5]). 
As the second moment a = (s 2 P(s))^ of the NNSD de- 
creases with increasing disorder, one could alternatively 
obtain an estimate for the roughness of a flake from a. 
We find however, that the dependence on (3 is numeri- 
cally more reliable, as the entire distribution is used for 
a fit to Pjj- 

While the shape of the unfolded NNSD, P(s), is sen- 
sitively dependent on disorder, it is to a good degree of 
approximation size independent, as our numerical data 
confirms (not shown) . The reason for this is the normal- 
ization (sP(s)) — 1, which scales out size-related effects. 



By contrast, the absolute level spacing (AE) contains di- 
rect information on the size of quantum confinement. By 
rescaling each level spacing by the local energy, relative 
to a fixed energy E$ = lOOmeV, one obtains the energy 
independent expectation value 

^ s2 2 5500nm 2 
= (^f) 2 — = d 2 rneV. (3) 
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This rescaled mean level spacing is, indeed, independent 
of edge roughness and disorder [see Fig. EJb)]. Agree- 
ment with the experimental data [|| for d = 40 nm is 
surprisingly good. 

In conclusion, the spectrum of realistic graphene quan- 
tum dots in the presence of disorder (edge roughness or 
defects) reveal unique features which differ from both 
Schrodinger or Dirac billiards of confined massive or 
massless free particles. The graphene bandstructure near 
the K point leaves clear imprints. They include inter- 
ference structures in the wavefunctions, enhanced con- 
finement effects, and a delayed transition from Poisson 
to Wigner-Dyson nearest-neighbor distributions. While 
one still "cannot hear the imperfect shape of the drum" , 
the size and roughness of graphene quantum dots can be, 
indeed, inferred from the spectral properties. 
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